Colloidal attraction induced by a temperature gradient 
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Colloidal crystals are of extreme importance for applied research, such as photonic crystals tech- 
nology, and for fundamental studies in statistical mechanics. Long range attractive interactions, 
such as capillary forces, can drive the spontaneous assembly of such mesoscopic ordered structures. 
However long range attractive forces are very rare in the colloidal realm. Here we report a novel 
strong and long ranged attraction induced by a thermal gradient in the presence of a wall. Switching 
on and off the thermal gradient we can rapidly and reversibly form stable hexagonal 2D crystals. We 
show that the observed attraction is hydrodynamic in nature and arises from thermal induced slip 
flow on particle surfaces. We used optical tweezers to directly measure the force law and compare 
it to an analytic prediction based on Stokes flow driven by Marangoni forces. 



Colloidal particles, dispersed in a homogeneous sol- 
vent, undergo unbiased Brownian motion unless an exter- 
nal body force, such as gravity, drives them along a par- 
ticular direction. Directed colloidal transport can also 
be induced, in the absence of a net body force, by a 
gradient in some external field such as electric potential 
(electrophoresis), temperature (thermophoresis) or con- 
centration of a solute molecule (diffusiophoresis) pQ. In 
phoretic transport, due to the absence of a net body force 
acting on the particle, and the thin interface layer around 
it, the surrounding fluid is perturbed with a flow field de- 
caying as 1/r 3 , compared with the leading term of order 
l/r,for sedimentation. As a consequence hydrodynamic 
interactions, which are crucial for sedimentation prob- 
lems, are usually negligible in phoretic transport. The 
opposite situation may occur when particles transport is 
impeded by an external wall. While sedimenting par- 
ticles become stationary as a result of a force balance 
between gravity and wall repulsion, phoretic particles 
have a zero velocity despite the presence of a non-zero 
body force (wall repulsion) acting on them. In such pe- 
culiar condition, gradients in the external field will pro- 
duce a flow field around a stationary particle which will 
be long ranged (1/r) determining strong hydrodynamic 
couplings to nearby particles. Such hydrodynamic in- 
teractions have been predicted by Squires [2 and are 
expected to manifest as an effective attractive potential 
when phoretic particles are pinned on a wall. Here we re- 
port the direct observation of strong attraction between 
colloidal particles confined at a solid interface by ther- 
mophoresis. The induced force is strong enough to lead 
to the formation of stable close packed structures and 
can be efficiently exploited to drive quick self assembly 
of colloidal crystal. Evidences for such a "hydrodynamic 
attraction" in temperature gradients have been indepen- 
dently reported by Weinert and Braun [3|. We propose 
an analytic expression for the effective force and validate 
it by direct force measurement with optical tweezers. 

A thin (~ 15 fim) sample cell is obtained by sand- 
wiching a small drop of sample solution (2 /am diameter 



silica beads in a glycerol/ water mixture) between a cov- 
erslip and an absorbing yellow filter as the microscope 
slide. High temperature gradients (~ \Kj\im) can be 
produced by focusing a laser beam on the absorbing ceil- 
ing of the sample cell (Fig. [I]). The same microscope 
objective (lOOx NA 1.4 in a Nikon TE2000U inverted mi- 
croscope) is used both for focusing the heating beam and 
for bright field imaging. As soon as the laser is turned 
on, the particles move towards colder regions under the 
effect of thermophoresis, following the direction of the 
negative thermal gradient, as expected for thermophobic 
behavior [4]. We used a finite element method (COM- 
SOL Multiphyics) to solve the stationary heat equation, 
with a heat source term modeled as a Gaussian beam, 
that propagates in an absorbing medium (top coverslide) . 
The distribution of the isothermal surfaces (black lines) 
and the temperature gradient direction (white lines), due 
to the absorption of a laser beam by the optical filter at 
the top, are depicted in Fig. 1. Once the particles reach 
the bottom wall, a strong interparticle interaction man- 
ifests leading to the formation of closed packed clusters. 
In order to reduce the radial components of temperature 
gradient we focused on the top absorbing wall an array 
of 20 traps along a circumference of 30 /am diameter. Iso- 
lated particles in the inner ring region are pushed against 
the bottom wall and then perform unbiased 2D Brownian 
motion. On the other hand thermophobic particles in the 
outer region are prevented to enter the heated region. In 
this way we can have a small number of particles in the 
inner region and record the dynamics of cluster forma- 
tion . Each cluster was able to find the closest packing 
structure in a few tens of seconds. A few sample clusters 
are shown in (Fig. [5J a). These clusters are completely 
reversible and disgregate as soon as the temperature gra- 
dient is turned off (Fig. [2j b). At higher concentrations 
particles assemble in stable, two dimensional crystal (Fig. 
3c). 

The existence of a confinement-induced long range at- 
traction between particles of like charge has been ob- 
served and debated during the past ten years [5j [6] . Such 
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FIG. 1: Temperature field in the cell. Distribution of the 
isothermal surfaces (black lines) and the temperature gradient 
direction (white lines), due to the absorption of a laser beam 
by the optical filter at the top. 




interaction is also found to increase when the like charged 
particles penetrate the electric double layer at the inter- 
face j5j [7]. One could think that the only role of the 
thermal gradient is that of pushing the particles against 
the double layer of the bottom interface and that the ob- 
served interactions are ultimately electrostatic in nature. 
However like-charge attraction disappears on a metallic 
surface [HE]. We then checked (Fig. [2] d) that the at- 
traction observed in our experiment persists on a metallic 
surface obtained through a gold coating 10 nm thick over 
the glass surface. The position of the laser beam inducing 
the thermal gradient is indicated by the brighter particle 
at the bottom, trapped by the beam and displaced from 
the bottom wall. Close-packed structures form both on 
a dielectric charged surface, i. e. the cover glass (visible 
as the clearer band), and on a metallic surface (darker 
band). Neither we can attribute the observed attraction 
to an imaging artifact [TO] , as the formation and stabil- 
ity of small clusters clearly evidences the presence of a 
strong attractive force. Another candidate for the origin 
of the attraction could be the distortion of the isothermal 
surfaces in the fluid, due to the different thermal con- 
ductivity of the particles and to the spatial asymmetry 
induced by presence of the wall. As, in our case, ther- 
mophoresis drives particles along the negative thermal 
gradient, this effect may induce an attraction between 
particles touching a colder wall, when thermal conduc- 
tivity is higher than the solvent's one, or a repulsion in 
the inverse case. However we always observe attraction 
in both cases of a higher (silica) and lower (polystyrene) 
thermal conductivities. 

We suggest, in accordance to [3], that the underlying 
mechanism is hydrodynamic in nature. The basic idea 
is that even if the particles stop when they reach the 
bottom wall there's a non zero net force acting on them, 
provided by electrostatic repulsion from the wall (the net 
force and torque produced by thermal gradients on the 
interface are always zero [lj). At large enough particle 
separations the ambient flow produced by one particle 
at the location of the other can be obtained as the flow 
propagating from a point force. In the absence of the wall 




FIG. 2: Reversible self-assembly of colloidal particles in a 
thermal gradient, a) Particles confined to the bottom wall by 
thermophoresis assemble in close-packed structures, b) Re- 
versible aggregation (first row, temperature gradient on) and 
disgregation (second row, temperature gradient off) of closed 
packed clusters, c) High concentrated samples assemble in 
stable, two dimensional crystal, d) Close-packed structures 
form both on a dielectric charged surface, i. e. the cover glass 
(visible as the clearer band), and on a metallic surface (darker 
band) obtained through a gold coating 10 nm thick over the 
glass surface. 



such a flow would be described by the so called Stokeslet 
[T2] which wouldn't have any attractive component on 
the other particle. When the no-slip condition imposed 
by the wall is taken into account, image singularities ap- 
pear (Stokeslet, source doublet and Stokes doublet [TT] ) 
contributing a flow component u x pulling the the two 
particles together. The streamlines of the flow produced 
by a Stokeslet pointing away from the wall are reported 
in Fig. [3] Such a flow will have, on the center of a second 
particle touching the wall at a distance s, a component 
parallel to the wall given by [2]: 
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FIG. 3: Streamlines of the flow generated by a point force 
pointing away from a plane wall with sticky boundary con- 
ditions. The resulting flow can be represented as the super- 
position of the Stokeslet flow produces by the force F plus 
an image singularity sister consisting of a Stokeslet, a Stokes 
dipole and a source dipole [11] . 



where F is the force exerted by the wall on the particle 
and h is the distance between the particle center and the 
wall. In the limit of large interparticle distances s (with 
respect to the distance h between particles' centers and 
the wall), the flow component parallel to the wall decays 
as 1/s 4 . The slowly decaying components (1/s 2 ) coming 
from the image Stokeslet and Stokes doublet cancel out 
at any s. Such an ambient flow on the second particle will 
be undistinguishable from an attractive force of intensity: 



f = 67TT]a\u x 
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where A = X(a/h) is an adimensional correction factor 
on the Stokes drag due to the presence of the wall [13]. 
The particles will then perform Langevin dynamics un- 
der the action of such interaction forces in addition to the 
usual drag and stochastic forces. Such an effective force 
field is derivable from a potential function and therefore 
particle statistics will be governed by Boltzmann statis- 
tics with an effective interaction potential: 
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The force F needed to hold a particle in a temperature 
gradient can be evaluated along the same lines suggested 
by Wurger [14] for the calculation of the thermophoretic 
speed U. In that paper the author solves Stokes equa- 
tions for a homogeneous fluid imposing the three bound- 
ary conditions given by: i) a vanishing net force is ap- 
plied to the particle, ii) the tangential component of the 
stresses on the particle surface cancel out the Marangoni 



force, iii) a zero normal component of fluid velocity rela- 
tive to the particle surface, which in the lab frame gives: 



U 



(4) 



In order to satisfy the above conditions one gets a value 
for the thermophoretic speed given by [T4] 



U = — VT 
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where r] is the solvent viscosity, = dj/dT is the par- 
tial derivative with respect to temperature of the interfa- 
cial energy 7, and n = 3k,s/(2k>s + ^p) is a dimensionless 
combination of particle and solvent thermal conductivi- 
ties. 

We seek here a solution for a spherical particle held 
stationary in a thermal gradient by a body force. The 
corresponding boundary conditions become: i) again a 
zero normal component of fluid flow on particle surface 
in the particle frame which now coincides with the lab 
frame, ii) a non vanishing, but unknown in modulus, net 
force F on the particle, iii) an off-diagonal stress ten- 
sor canceling the Marangoni force. A Stokeslet flow has 
a zero off diagonal stress tensor, a non vanishing body 
force F, and a fluid flow projection on the particle sur- 
face normal ft given by 



(6) 
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Comparing Q to ([6| it is easy to satisfy the above 
conditions by simply adding to the flow obtained above 
for the freely moving particle, a Stokeslet flow with 



F = -47T77U = — VT = -Vr (7) 
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where T is the total interfacial energy. We can now 
assume that F remains a good estimate even in the pres- 
ence of the wall and get an expression for the force / by 
substituting ^ in ([2|. For spheres almost touching the 
wall (h ~ a [15]) we get: 
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The force depends linearly on the interface energy gra- 
dient, increases with the square of particle size but is 
independent from solvent viscosity. 

In order to validate the above expression we used opti- 
cal micromanipulation to directly measure the force law 
with interparticle distance. 
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Holographic optical tweezers allow to isolate a single 
pair of particles and vary their relative distance while 
gauging their interaction. In particular, a high intensity 
optical trap (P = 42 mW) is focused next to the interface 
with the top wall in order to produce the thermal gra- 
dient, while two traps of lower intensities (P = 3 mW) 
are used to hold two particles on the bottom wall. We 
checked that the intensity of the last traps alone is low 
enough not to induce further attraction between the par- 
ticles. In order to avoid fluid instabilities induced by the 
high thermal gradient [16 , we have chosen a high viscos- 
ity suspension (a glycerol-water mixture 56% w) and a 
thin cell gap (13 fim thick). 

The thermal induced interaction tend to push the par- 
ticles towards each other and out of the optical traps, un- 
til the restoring trap force balances the attractive force. 
By measuring the difference As between trap centers dis- 
tance and interparticle distance, the value of the attrac- 
tive force can be deduced: / = kAs, where k is the equiv- 
alent elastic strength of the two traps system p~7 and 
may be determined experimentally from the mean square 
displacement of particle distance: (Ad 2 ) = KsT/k. In 
the reported experimental conditions we found k — 0.1 
pN//im. However an experimental determination of the 
prefactor f in ^ is very hard because of the strong 
dependence of A on the exact value of h and difficulties 
in measuring the temperature gradient. On the other 
side the dependence on interparticle distance only con- 
tains particle radius and can be accurately checked. Fig. 
[4] reports the experimental determination of force as a 
function of interparticle distance together with a fit to 
Eq. [8] with /q = 4.8 pN as the only fit parameter. 

In conclusion, we have shown that a thermal gradi- 
ent pushing thermophoretic particles against an inter- 
face can induce a strong and long ranged attraction be- 
tween like-charged particles, leading to the prompt for- 
mation of stable, ordered structures. We provide a static 
measurement of this interaction and quantitatively con- 
firm the prediction of a hydrodynamic pseudo-force [2] 
driven by Marangoni forces on particle surface. The long 
range nature of the interaction, similarly to capillary 
forces, may open new routes to fabrication of ordered 
structures. Already, this mechanism may have enhanced 
the self-assembling observed under a convective flow and 
uniquely attributed to a confining effect [18, 19| [20]. 

We thank M. Ortolani and the IPN-CNR labs for 
kindly providing us the gold coated cover glass. 
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FIG. 4: Intensity of the attractive force as a function of 
distance. Open symbols are experimental data, solid line is a 
fit to Eq. J8l. 
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